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ABSTRACT 

Several classes of stars (most notably O and B main-sequence stars, as well as accreting white dwarfs 
and neutron stars) rotate quite rapidly, at spin frequencies greater than the typical g-mode frequencies. 
We discuss how rapid rotation modifies the K-mechanism excitation and observability of g-mode oscilla- 
tions. We find that, by affecting the timescale match between the mode period and the thermal time at 
the driving zone, rapid rotation stabilizes some of the g-modes that are excited in a non-rotating star, 
and, conversely, excites g-modes that are damped in absence of rotation. The fluid velocities and tem- 
perature perturbations are strongly concentrated near the equator for most g-modes in rapidly rotating 
stars, which means that a favorable viewing angle may be required to observe the pulsations. Moreover, 
the stability of modes of the same I but different m is affected differently by rotation. We illustrate this 
by considering g-modes in Slowly Pulsating B-type stars as a function of the rotation rate. 

Subject headings: stars: early-type — stars : oscillations — stars : rotation — stars: variables: other 
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1. INTRODUCTION 



Despite a long-standing interpretation of multiple peri- 
odicity in Ce pheids and Slowly Pulsating B-type (SPB, 



Waelkens 1991 ) stars (also r eferred to as 53 Per st ars, after 
the prototype discovered by |Smith fc Karp 1976|) in terms 
of nonradial pulsations, the details of the excitation mech- 
anism remained unresolved until quite recently. The driv- 
ing of pulsations in the He + ionization zone, which gives 
rise to the classical Cepheid instability strip, cannot oper- 
ate in these hot stars. In a remarkable example of the in- 
terplay between t heory and observation, the updat e d opac- 
ity calculatio ns ( Iglesias, Rogers, fc Wilson 1992| ; Bcatorj 



et al. 1994 ) found a significant opacity increase in the 
metal partial ionization zone, just as suggested by Simori| 



1982) as the explanation for intrinsic variability in 53 Per 
stars. The resulting instability region, determined by the 
k mechanism operating in the "metal opacity bump" at 
T « 1.5 x 10 5 K, incorporates most OB stars (M £ 3Af©) 
that are in the core hy drogen burning stage ( Dziembowskj 
fc Pamyatnykh 1993 ; Dzicmbowski, Moskalik, & Pamy- 
atnykh 1993; Gautschy & Saio 1993). The most recent 



excitation calculations for Cepheids and SPBs ( Pamyat 



nykh 199S ) place all known variable stars in these classes 
within the theoretical instability region. However, it is not 
clear that all stars in this region are indeed pulsating. 

All previous g-mode pulsation calculations have been 
carried out for non-rotating or very slowly rotating stars. 
However, most OB stars rotate relatively rapidly and ob- 
servations indicate that g-mode pulsations and rotation 
are somehow related. As discussed in §||, rotation strongly 
affects g-mode properties when 2P mo( jo > P S pin, where 
-Pmode is the mode period in the co-rotating frame, and 
P sp in is the spin period. Rotation would have to be near 
breakup to affect the p-mode pulsations in the (3 Cepheids. 
However, typical SPB stars (M w 4M©, R xt 3i2 Q ) exhibit 
high radial order g-mode pulsations of day-like periods, so 



that rotational velocities 



v ro t > 75 km s 



1 day 

-fmodc 



R 

3R7. 



(1) 



will significantly affect the g-modes. The mean rotat ional 
velocity for B stars is v rot ~ 200 km s _1 (Lang 1991) and 
thus cannot be neglected in pulsation calculations. 

One of the hallmarks of the photometrically discovered 
SPBs is that they are anomalously slow rotators, with 
projected rotational velocities wsini ^5 — 40 km s _1 

70 



(Waelkens 1987; 


Waelkens 1991 


. While most of the r 


SPBs discovered by Hipparcos ( 


Waelkens et al. 1998) 



also slow rotators, some have usini's up to 150 km s _1 , 
a regime where rotational effects are very important. An- 
other correlation between rotation and pulsations is that 
no SPB stars were found in the open clusters NGC 3293 
and NGC 4755, where the mean projecte d velocity of B 
stars is wsini w 170 km s _1 ( |Balona 1994a ). Nevertheless, 



there were many Cep stars discovered •photometrically i n 



the same clusters (Balona 1994b; Balona & Koen 1994), 



so this was not due to a selection or mctallicity effect. The 
fact that most photometrically discovered SPBs are slow 
rotators, as well as the a bsence of SPBs in these two clus- 
ters, led Balona (1995)| to conjecture that rapid rotation 



suppresses g-modes in these stars. 

In this paper we explore the interaction between rapid 
rotation and nonradial pulsations. We begin in by out- 
lining the effects of rotation on the adiabatic mode proper- 
ties (frequencies and eigenfunctions) and observability of 
g-modes in rotating B stars. In §[| we discuss the interplay 
between rapid rotation and excitation of pulsations by the 
k mechanism and its effect on the stability of g-modes in 
B stars. We close by outlining the future directions for 
our work and mention other potential applications of the 
rotationally modified g-mode theory we are developing. 



2. ADIABATIC G-MODES IN RAPIDLY ROTATING STARS 



1 



2 



Rapid Rotation and g-Modes in B stars 



We begin by reviewing the theo ry of adiabatic nonra- 



approximation, 



dial p ulsations in rotating stars (see Bildstcn, Ushomirsky^ 
& Cu |bler 1996| , hereafter BUC, for a detailed discussion 
and references). Throughout this paper, we assume that 
the star rotates uniformly with spin frequency fi, much 
smaller than the breakup frequency fib « (GAf/i? 3 ) 1 / 2 , 



and ai lway3 work in the co - rotating frame. — The atar is 



nearl y ephorical in this limit, and the offoctc of the con — 

trifugal force (which are of 0(fl/flb) 2 ) can be neglected. 
For our purposes, "rapid rotation" is the regime in which 
uj < fl <C fib, where w is the mode frequency in the co- 
rotating frame. In this case, the primary difference in 
the momentum equation is the Coriolis force, the effect 
of which is of 0(fl/u>), and is always fully included in our 
calculations. Perturbation expansions typically employed 
in the slow rotation regime break down when rotation is 
rapid. We use the method developed by BUC which con- 
verges regardless of the rotation rate. When fl ;> u>, the 
Coriolis force strongly affects the g-modes. On the other 
hand, neither the p-modes nor the radial modes (both of 



which have w > ft),) arc strongly modified when fl fib- — 
Wic hin the traditional and — Cowling approximations 

(see BUC and references therein), the adiabatic pulsation 
equations separate, i.e. the solutions can be written as 
Q{r,9,<j),t) = Q r (r)Q e (e)e im < l ' +luJt . The frequency in the 
observer's frame is then luj = u> — mfl. The Cowling ap- 
proximation (i.e. neglecting the perturbation of the grav- 
itational potential) is justified because we are studying 
high radial order g-modes. The traditional approxima- 
tion amounts to neglecting the components of the Coriolis 
force due to the local horizontal component of the spin. 
For g-modes, this is justified when the scale height h is 
small compared to the radius, h/r <C 1, and the buoyancy 
is strong compared to the radial component of the Cori- 
olis force, |iV 2 | 3> rfloj/h, where N is the Brunt- Vaisala 
frequency. Both of these conditions are satisfied in the 
envelopes of B stars and allow us to make some progress 
without resorting to a 2-D calculation. Extending our cal- 
culation into the deep interiors will force us to drop this 
approximation, however, since most excitation and damp- 
ing is confined to the envelope, the results presented in 
§0 are qualitatively correct. In the traditional approxi- 
mation, the angular part of the solution obeys Laplace's 
tidal equation, L(Qg) = —XQe, where L is an angular 
operator that depends only on m and the rotation param- 
eter q = 2fl/uj (BUC). Typical solutions X(q) for g-modes 
are shown in Figure 2 of BUC. The radial equations are 
identical to those for a non-rotating star, except with the 
spherical harmonic eigenvalue 1(1 + 1) replaced by A. When 
buoyancy is strong, ./V 2 3> uj 2 , if luq is the frequency of an 
I = 1 mode in a non-rotating star, then the frequency of 
the same mode at arbitrary spin is just uj — luq(X/2) 1 / 2 . 
Moreover, when TV 2 » uj 2 , the radial dependence of the 
eigenfunctions is the same is in the non-rotating case. 

One can view A 1 / 2 as the effective transverse wave num- 
ber (i.e. fc t 2 , = X/R 2 ). For no rotation (q = 0), A = 1(1 + 1) 
and is independent of m. For slow rotation (g < 1), the 
degeneracy between the different m values is lifted, lead- 
ing to the familiar m-splitting. However, for rapid rota- 
tion (q ^> 1), the modes with > 1 nodes in the angular 
eigenfunction have A ~ (2/ M — l) 2 q 2 . Hence, when the star 
is rapidly rotating, the eigenfrequency in the traditional 



,1/2 



(2Z M - l)fiw , 



(2) 



is roughly independent of m, and is just a geometric mean 
of the spin frequency fl and the no n-rotating I ~ 1 mode 



frequen cy loq for modes with L > 1. Papaloizou & Pringle 



1978)| first found this scaling in the WKB limit (Z M > 1 



and BUC extended their result to the observationally rele- 
vant low-? modes. On the other hand, X(q) for modes with 
1^ = zero crossings in the angular direction (i.e. the 
modes with m = —I) stays low and approaches A ~ m 2 for 
q > 1. 

The angular dependence of the eigenfunction for rapid 
rotation is no longer just a Yi m , however, we label the 
g-mode at arbitrary fl by the I and m values the same g- 
mode has when fl — 0. The fluid displacements and pres- 
sure perturbations for the > 1 modes are strongly con- 
centrated towards the equator, with the bulk of the ampli- 
tude bet ween the latitudes —l/q< cos 9 < 1/q (Figure jj]; 
see also |Lee & Saio 1990|; BUC; Figure 5 of (Townsenc] 



1997). Therefore, g-modes in rapidly rotating stars may 



be harder to detect because pulsation amplitudes are ap- 
preciable only over a small fraction of the stellar surface. 
How this affects detectability depends on how large the 
amplitudes are for the modes excited in the rapidly rotat- 
ing case compared to the non-rotating case. 

3. EXCITATION OF ROTATIONALLY MODIFIED G-MODES 

The most prominent effect of rapid rotation (5 > 1) on 
g-mode pulsations is that their periods in the co-rotating 
frame become shorter, -P mo dc oc (-Pspin-fo) 1 ^ 2 , where Pnis 



the period of the mode in a non-rotating star (see eq. |2| ) . 
This scaling is crucial for understanding the effect of rapid 
rotation on mode stability, as we now show. 

The qualitative understanding of the effect of rotation 
on g-mode excitation comes from the work integral formu- 
lation. Aizcnman & Cox (1975) derived the work integral 
for a star possessing arbitrary fluid flow in the equilib- 
rium configuration. Since we assumed rigid rotation, the 
background stellar model cannot be in thermal balance. 
However, the thermal imbalance effects are of 0(fl/flb) 2 
and are negligible compared to the dominant K-mechanism 
terms (Aizenman & Cox 1975). The dominant contribu- 



tion to the work integral is the familiar expression, 



A ( , -F\dM = J ^dr, (3) 



where AT is the Lagra ngian temperature perturbation 
(Gautschy & Saio 1995), and the second equality is ob- 
tained by integrating over the angles. A mode is lin- 
early unstable if W > 0, so that regions of the star with 
dW/dr > excite the pulsation. Driving occurs in re- 



gions with d[ttT 



((91og«/G>lo gr) and k 



>/(r 3 



l)]/dr > 0, wh ere kt = 
p — (d log n/d log p)t (Gautschy 
fc Saio 1995| ) , and competes with damping throughout the 
rest of the star, where the above inequality is reversed. 

We consider pulsational stability of several g-modes in 
a 4Mq main-sequence model (t »s 10 8 years, approxi- 
mately half of main-sequence lifetime) with log L/Lq — 
2.51, logTegf = 4.13, and Z — 0.02, computed using the 
Warsaw-New Jersey stellar evolution code provided to us 
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by R. Si enkiewicz. 



used by Dziembowski ct aL (1993) 



This model is ver y similar to the one 
We use the quasi- 



adiabatic approximation, i.e. evaluate W using the adi- 
abatic eigenfunctions. At some point in the stellar en- 
velope the local thermal time, T th (r) = J^r p TrtM/T, sa 



CpTA\M-/-L, becomes much smaller than the mode period 



Prr 



These nonadiabatic parts of the envelope do not 



contribute to either driving or damping, so we truncate the 
work integral at Tth = Pmodc to approximate this effect. 

In SPB stars, the k mechanism acts in the metal ioniza- 
tion zone, at logT i=s 5.2. Figure ^a shows Kt + K p /(T3 — 1) 
for the AMq SPB model described above in the region 
of the metal opacity bump. Figures ||b and ||c show the 
differential contribution to the work integral, dW/d\ogT, 
for I = 2, n = 14 and I = 1, n = 30 modes, evaluated 
using adiabatic eigenfunctions. The contribution to the 
work integral from the inner part of the ionization zone 
is positive, while the contribution from the outer part is 
negative, a typical behavior for the g-modes in SPB stars. 
As n increases, P mo de in the co-rotating frame becomes 
shorter. However, as discussed in §||, the radial depen- 
dence of the eigenfunctions in the traditional approxima- 
tion is the same as in the non-rotating star. Thus, for 
\ow-l SPB pulsations, the dominant effect of rotation is 
to truncate the work integral closer to the stellar surface. 
For the mode in Figure |2|b, the period in a non-rotating 
star is P mo dc = 0.675 days, so the work integral is trun- 
cated at logT sa 5.18, resulting in an unstable mode, 
W > (see Figure ||). However, as the spin frequency 
is increased from O = 0, the period of the mode decreases, 
and the location in the star where the mode goes nona- 



diabatic, Pry 



}dc 



r t h, moves closer to the stellar surface. 



In the non-rotating star, all of the ionization zone already 
contributes to the excitation of the mode, so at moder- 
ate rotation, the outer, damping, part of the ionization 
zone starts contributing as well. As shown in Figure [| 
at P S pin ~ 30 hours (w ro t ~ 140 km s _1 ), the m = 1 and 
m = 2 harmonics of this mode become stable. On the 
other hand, the mode in Figure ||c, with a 2.72 day pe- 
riod is stable when SI = 0. However, as rotation becomes 
rapid, the adiabatic part of the star encompasses more of 
the driving part of the ionization zone. As shown in Fig- 
ure H, the m — 1 harmonic becomes unstable (W > 0) 
when P S pi n ~ 70 hours (t> ro t ~ 60 km s _1 ). 

As evident from the preceding discussion, if the ion- 
ization zone is completely adiabatic (i.e. if the thermal 
time at the ionization zone Ti on iz 3 s Pmode), the driving 
due to the inner part of the ionization zone is likely to be 
cancelled by the damping due to the outer part. On the 
other hand, if the ionization zone is completely nonadia- 
batic (i.e. if Ti on i Z <C Pmode), then it does not contribute to 
either driving or damping. Thus, to get the optimal driv- 
ing contribution from the metal opacity bump in SPBs, 
the period of the mode in the co-rotating frame should be 
comparable to t he thermal time in the ionization zone, i.e. 



P mo do ( pautschy fc Saio 1993|) 



In close ag reement with the calculations of Dzicmbowskj 



et al 
and n 



(1993; 



we find that the 



14-23 modes with I = 1 



14-29, 1 — 2 modes are excited in the non-rotating 



4M Q SPB model. Rapid rotation damps the modes with 
I = 1, n = 14-16 and I = 2, n = 14-16, and excites the 
I = 1, n = 24-39 and I = 2, n = 30-40 modes. On the 



other hand, for this model, the modes with n ~ 20, in the 
middle of the excited range, can be stabilized only if rota- 
tion is near breakup, or cannot be stabilized at all. This 
behavior can be understood in terms of the normalized 



growth ra te of the mode, r\ — W/ J \dW/dr\dr (Sterling- 



werf 1978), which measures the robustness of mode exci- 
tation. The modes with n rs 20 have r\ « 0.1-0.2, while 
the ones with n w 15 or 30, at the boundaries of the ex- 
cited range, have ij « 0.01-0.05 (Dziembowski et al. 1993). 
This means that if rotation is to stabilize the modes with 
n ~ 20, it has to change W by about 15-20%, which is 
difficult to achieve within our current approximations. On 
the other hand, for n m 15 or 30, a few percent change 
in W can change the character of the mode, and is easily 
accomplished by rotation. 

The modes with m > are more strongly affected by 
rotation than the m < modes, because the former reach 
larger A (shorter period) at smaller q (slower rotation). 
This is an interesting correlation that may be confirmed 
observationally. The I = 1, m = — 1 and I = 2, m = — 2 
modes (not shown in Figure |3|) are not strongly affected 
by rotation. The value of A for these modes only changes 
from 2 to 1 and 6 to 4, respectively, as rotation becomes 
rapid, and thus the periods and stability of these modes 
are not significantly modified. 

4. CONCLUSIONS 

Though we have made progress in understanding the 
interaction between rotation and pulsations, we have yet 
to unambiguously explain the preponderance of slow ro- 
tators among the field SPB stars. Rotational modifica- 
tion of the mode period can destroy the timescale match 
Pmodo ~ tioniz, and stabilize a g-mode that is excited in 
a non-rotating SPB model. Conversely, if the thermal 
time at the driving zone T; on i z does not match the non- 
rotating mode period, the mode can be "spun up" and 
become excited if the star is rotating rapidly. The periods 
of excited modes in the co-rotating frame are largely de- 
termined by the structure of the star, in order to satisfy 



Pn 



)de 



However, the radial orders n of the excited 



g-modes depend on rotation. In order to maintain the 
timescale match with increasing J7, one needs to increase 
n (note that the co-rotating frame periods and stability 
of g-modes with m — — I is not significantly affected by 
rotation). However, going to a larger n may mean more 
radiative damping in the stellar interior, especially in the 
/x-gradient zone. In addition, we find that different m har- 
monics of the same mode are affected to various degrees 
by rotation, which may have consequences for line profile 
observations capable of determining the I and m values of 
the pulsation modes. 

Rotation does not significantly affect the g-modes in B 
stars until the spin period is comparable to the mode pe- 
riod. Rotation this rapid confines most g-modes to the 
equatorial region, requiring a favorable viewing geometry 
for detection, and can destroy the timescale match needed 
for mode excitation. In this sense, the correlation between 
slow rotation and SPB pulsations may be an observational 
selection effect after all, since the g-modes in rapid rota- 
tors, even if present, only cover a small fraction of the 
stellar surface. In order to differentiate between these two 
effects, one needs information about the inclination an- 
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gle i, as well as the amplitude of pulsations. Though, as 
explained in rotation is ne ver "rapid" (in the sense 



2S1 ;> lu) for p-mode pulsations, Lee fc Baraffe (1995) have 
shown cases where the low-order rotational perturbations 
can affect p-mode stability in f3 Cep-like stars rotating 
close to breakup. 

Our work may also be applicable to other classes of vari- 
able stars, most notably the 7 Doradus st ars. These stars. 



name d afte r the prototype discovered by Cousins & War 
ren (l |963)| , are early F-type main-sequence stars that lie 
at or below the cool edge of the classical Cepheid insta- 
bility strirj J _and_show variability with periods of around a 
day (Breger 1995). It is currently believed th at the vari- 
ability is due to g- modes of high radial order ( KrisciunaE 



1997). For an F-type star with R ~ 1.5i?o, rotational ve- 



locities that would affect g-modes with one-day period are 
v TO t > 40 km s _1 (see eq. Q). Typical rotat ional veloc i- 
ties for early F stars are v TO t « 100 km s" 1 (Lang 1991). 



Therefore, if g-modes are indeed responsible for 7 Dor vari- 
ability, rotation will have a significant effect on the mode 
properties and stability. 
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Fig. 1. — (a) Adiabatic Eulerian pressure perturbations (colored contours) and fluid displacements (black arrows) for an 
I = 2, m = —1 g-mode in a non-rotating star. Normalization is arbitrary, but orange/yellow (purple/blue) contours denote 
positive (negative) pressure perturbations, (b) Same as (a), but for an I = 2, m = — 1 g-mode in a rapidly rotating star, 
with a spin period P sp i n = 0.5P mo d c - Note how pressure perturbations and fluid displacements are strongly concentrated 
towards the equator. 
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Fig. 2. — (a) The opacity derivative kt + K p /(T3 — 1) in a 4M Q SPB model, (b) Differential work integral, dW/d\ogT, 
evaluated using adiabatic eigenf unctions, for an I = 2, n = 14 mode. Dashed and dotted vertical lines show where r t h is 
equal to P m odc = 0.675 days, corresponding to the I — 2 mode period in the non-rotating star, and the breakup period, 
Pb ~ 8.6 hours, respectively. In the non-rotating case, dW/dlogT computed using nonadiabatic eigenf unctions would 
decay to zero for temperatures lower than that indicated by the dashed lines. We approximate this effect by truncating 
the work integral as described in §||. (c) Same as (b), but for an I = 1, n = 30 mode with P mo dc = 2.72 days for I = 1 in 
a non-rotating star. 
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Fig. 3. — Work integral W (eq. ||) as a function of the spin period P S pin (bottom axis) or equatorial rotational velocity 
w ro t (top axis) for our SPB model. Lines asymptoting to positive W for no rotation correspond to the I = 2, n = 14 mode 
shown in Figure while the lines asymptoting to W < for no rotation correspond to the I = 1, n = 30 mode shown 
in Figure ^p. Different m values are marked on the plot, where m > modes are retrograde. 



